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Abstract. We apply systematic methods previously used by Mihalkovic et at 
to predict the structure of the 'basic' Co-rich modification of the decagonal 
Al7()Co2oNiio layered quasicrystal, based on known lattice constants and 
previously calculated pair potentials. The modelling is based on Penrose tile 
decoration and uses Monte Carlo annealing to discover the dominant motifs, which 
are converted into rules for another level of description. The result is a network 
of edge-sharing large decagons on a binary tiling of edge 10.5 A. A detailed 
analysis is given of the instability of a four-layer structure towards c-doubling 
and puckering of the atoms out of the layers, which is applied to explain the 
(pentagonal) orientational order. 



1. Introduction 

A fundamental problem of theoretical solid-state chemistry is, assuming the ability 
to compute perfectly the total energy of a complex compound, to discover which 
arrangement of the atoms minimizes that energy. Quasicrystals offer a wonderful 
opportunity to face this problem, as many different arrangements are low in energy, 
but this approximate degeneracy is resolved by more subtle interactions, often in 
several steps of a hierarchy of length scales. 

Consider decagonal Al-Co-Ni, the pre-eminent decagonal Al-TM quasicrystal 
(TM= transition metal), a phase which is divided into many modifications of 
which the simplest are 'basic' Ni-rich (around Al7oCoinNi2o) and 'basic' Co-rich 
(around AlyoCo2oNiio). Mihalkovic et al [2 1 proposed a general method to discover 
the ground-state structure, taking as inputs the interatomic pair potentials, derived 
by Moriarty and Widom |3| using 'generalized pseudopotential theory', and the 
experimentally determined lattice constants. One initially describes the quasicrystal 
by a Penrose tiling of edge ao ~ 2.45 A with properly placed candidate sites on the 
tiles. Monte Carlo annealing, in which atoms hop on these sites, is used to find low- 
energy configurations. (This approach was used earlier by Cockayne and Widom [4 1 
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for decagonal Al-Co; Mihalkovic et al |2| added tile rearrangements as a helpful 
Monte Carlo move.)- In the next stage, motifs abstracted from these results are 
converted into rules to decorate larger tiles, thereby reducing the degrees of freedom 
and permitting larger simulation cells in which further orderings can be seen. 

2. Structure motifs and idealized decoration 

We have applied this method to the 'basic' Co-rich composition (more complete 
accounts may be found in Refs. From the initial stage simulation it became 

clear that the dominant motif is a decagon with edge ra Q « 4A and hence diameter 
2r 2 ao ~ 12. 8A; in this rest of this paper, the term ' 13 A Decagon' refers exclusively 
to this object. The exact occupancies (fig. ^) were clarified in the second stage 
of simulations. These were based on rao-edge rhombi in arrangments favoured the 
densest possible packing of non-overlapping 13A decagons made from the rhombi. 
(Overlaps were verified to be unfavourable energetically.) The atoms in a cluster 
are grouped into concentric rings (see figure^). Just inside the perimeter, a few 
Al sites form a partial 'ring 2.5' with elusively context-dependent occupancies. 
(The Al positions on each 13A Decagon edge (ring 3) also may depend on the 
13A Decagon's environment.) Almost all the TM atoms in the 13A Decagon are 
Co, which is presumably why it emerges at higher Co content. (This was also 
discovered by Hiramatsu and Ishii |7|, who did an exploratory study using the same 
code of a wider range of AlCoNi compositions.) 

The best arrangements for edge sharing of the 13A Decagons amounts to a 
10.5 A edge 'binary tiling' |8 9|, with the 13A Decagons sitting on the Large- 
atom sites of the binary tiling. Filling the space between is a secondary motif, which 
is centred on the Small-atom sites, which we call the 'Star cluster' and show in 
figure Qb). The centre of this cluster has a pentagon of five candidate TM sites, 
which are variably occupied by TM (almost always Ni) or Al. 

The TM occupancies show context dependence, which we will approximate 
with a very simple rule (figure^). We have four classes of TM atom. On each of the 
binary tiling's Large vertices we have a 13A Decagon with five Co(l) in its centre 
and ten TM(3) on its perimeter (ring 3). (For visual clarity, the TM(3) are not shown 
in figure^.) Now, on each (binary-tiling) edge connecting a 13A Decagon to a Star 
cluster centre, a Star cluster TM is strongly favoured or disfavoured according to 
whether the 13A Decagon's first ring has Al or TM on the same edge. (The distance 
between these sites is b 1 ' 2 ^ w 5 .6A, which is a maximum of Vtm-tmM |3|.] So, 
on every binary-tiling rhombus exactly half the edges do not pass over a Co(l): let 
us place on each such edge a Ni(s) within the Star cluster, and make the TM(3) site 
a Ni(3); all other TM(3) become Co(3). (As in the 'basic' Ni-rich composition, the 
difference between Ni and Co potentials favours Ni occupancy in TM-TM pairs and 
Co in lone TM sites.) 

Finally, as this too is energetically favoured, we place 2 Ni(s') wherever two 
Star clusters overlap, that is once on each Thin rhombus of the binary tiling. In 
a decagonal tiling, there is an average of 3.19 Ni on each Star cluster and the 
composition is AZ7oCo2i.2Nig.8 (assuming the TM fraction is 30%). (One other 
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Figure 1. Decoration model (a). Decoration of a 13A Decagon, named in 
concentric rings. Size of circles indicates which of the two layers; some third- 
ring 'Co' atoms become Ni (as shown) but not many. Dashed ovals mark four (out 
of eight) of the places showing Co-Al-Al-Co in projection, which after relaxation 
become a channel like figure |2^. (b). Decoration of Star cluster. There are 
additional Al along the outer edges, (c). 13A Decagons and Star clusters shown 
with dashed edges, placed on vertices of the binary-tiling rhombi (edge T 2 ao ~ 
10. 5A). Filled circles are TM atoms in 13A Decagon and Star cluster interior. 
The ± sign on each 13A Decagonindicates Ui = ±1. Ni(3) sites are shown in 
grey; there are Co atoms (not shown) on all other corners of the 13A Decagons. 



version of the idealization, including the placement of Al atoms, is described by 
Gu et al 1 6 1, and gave a composition Al70.1C022.4Ni7. 5, and it is even easier to devise 
variants with higher Ni content.) 

This model does not include the 'pentagonal bipyramid' motif that is prominent 
in the related approximant W( Al-Co-Ni) 1 10 1, so we have not addressed whether that 
motif is favoured in the decagonal Al-Co-Ni structure II II . 

3. Orientational order and puckering 

A full structure description must specify, in each 13A Decagon's first ring, the 
orientation of the pentagon formed by its ring 1 TM atoms (this is prominent in 
electron microscope images). Each orientation is parametrized by a spin-like variable 
a = ±1, with (7 = +1 when the C05 pentagon (projected on the xy plane) points 
up, that iswhen the five Co lie in the even-numbered atom layers. We seek an 
effective interaction of neighboring 13A Decagon's, like the spin interaction in an 
Ising model. A cluster arrangement with all clusters oriented the same will be called 
'ferromagnetic'; that quasicrystal has global pentagonal (not decagonal) symmetry, 
as is observed in fact for one of the Co-rich modifications 1 121 1131 . 

3.1. Fixed site list: interaction via Star clusters 

When we simulate with a fixed site list, there are small energy differences that favour 
either the 'ferromagnetic' or 'antiferromagnetic' order of 13A Decagon orientations, 
depending on the atom density assumed and which approximant is studied, showing 
there are competing contributions to the cluster-cluster effective interaction. (The 



Co-rich decagonal Al-Co-Ni 



4 



'ferromagnetic' effect wins for 0.068 — 0.074 A~ 3 , which includes the physically 
more reasonable densities.) 

Where does the interaction come from? The only atoms that change with a 
cluster's orientation are the central Al and the ring 1 AI5C05, but those are too distant 
to interact directly. (Our potentials were cut off around 7 A.) An indirect interaction 
may be mediated by other atoms in two ways. Firstly, the variable Al atoms in 
rings 2.5 and 3 are within range of ring 1 of both clusters; this guides whether the 
ring 3 Al (in projection) divide the 13A Decagon edge in the ratio r _1 : r~ 2 or 
t~ 2 : r _1 (See figure^.) This contribution favours 'ferromagnetic' order, so the 
shared 13A Decagon edges will match. (Expressed in energy terms, the argument 
is related to the one for the 'puckering' case in Sec.0] below.) 

Alternatively, neighboring TM atoms in the Star clusters interact with each 
other and also respond to the first-ring orientations of adjoining 13A Decagons: 
we suggest this is the origin of the 'antiferromagnetic' coupling, based on the TM 
rules of Sec. [2] Let us assume a repulsive energy V72 for Ni-Ni pairs on edges 
differing by 72° as found by Mihalkovic et al Q in 'basic' Ni modification. One 
such pair appears on each Fat rhombus having 13A Decagons with the same 
orientation. The effective orientation interaction is then s(l + 01 o^)!^ favouring 
the 'antiferromagnetic' arrangement. Thus, the two indirect effects are indeed 
competing. 

4. Relaxation and puckering 

Relaxation and molecular dynamics of the structures found in the fixed-site 
simulation were used to explore subtler aspects corresponding to small energy 
differences - including the relative orientation of neighboring clusters, which is a 
major focus of this paper. The most striking effect of letting atoms off of the fixed 
sites was that Al atoms in ring 2.5/ring 3 tended to run to new locations, in which they 
displaced from the layers in the z direction (puckering of atom layer). When we used 
a simulation cell with four layers, pattern of puckering is such that lattice constant 
doubled to 2c w 8 A. (That is the observed period of most Co-rich decagonal Al-Co- 
Ni modifications and approximants, and it was long known crystallographically that 
puckerings are the major way in which the real cells are doubled.) 

4.1. Channels 

How shall we best understand the positions taken by Al after they are freed from 
the fixed site list? Let us consider the net potential felt by an Al given the TM 
atoms, which we fix on ideal sites (they hardly displace under relaxation in any case). 
Winding one-dimensional 'channels' are found along which the potential is low and 
comparatively fiat, typically between two columns of Co atoms (a ubiquitous pattern 
in the structure) forming a zigzag ladder, as shown in figure^. (Henley et al | 14 1 
discovered similar channels in 'basic' Ni-rich modifications from the time-averaged 
Al density in a molecular dynamics simulation.) 
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Figure 2. 'Channels' accomodating three Al atoms per four layers, (a). Path of a 
channel favourable to Al between two staggered columns of Co atoms; positions of 
the rigid site list shown as dashed circles. The farther TMs (grey) actually lie in a 
direction nearly normal to the paper, (b). Single-Al potential U (z) in the channel, 
showing minimum-energy configuration of three atoms with period 2c. (c). Four 
channels adjacent to single Co column; the two channels in the plane normal to the 
paper break the period c symmetry in the opposite way from the other two. 



The single-Al potential in a channel may be approximated as 

~ /47rz\ ~ / 2irz \ 

U(z) =U - [/ c/2 cos(— ) +a£/ c cos(— ) (1) 

Using our potentials, U c ~ U c /2 ~ 0.2 eV is found. Let us explain the origin of these 
coefficients. The Co atoms adjacent to the channel have symmetry under one-layer 
shifts, which explains the period-c/2 term of Q. Its minimum occurs in each atom 
layer, where the Al atom would be equidistant from three Co atoms at r ~ 2. 5 A, the 
minimum of the VaiCo (t) pair potential. 

Figure [2] shows another column of TM farther away which breaks the 
degeneracy between even and odd layers. This explains the period-c term in equation 
0, in which u is the 'spin' variable defined above, and a = +1 (resp. —1) when 
the distant Co atoms are in even (or odd respectively) layers. (The deeper minima 
occur in the TM-poor layer.) Note that U c is proportional to — oJVaitm (r) /dr, which 
happens to have a maximum at r « 4A, the distance to the distant column. 

4.2. Puckering 

Now, three Al atoms can and do fit into every four layers: their mean separation 
2c/3 k 2.7k is right at the Al-Al hardcore radius. Just where do they prefer to lie in 
each channel? The optimum configuration has a symmetry such that z% m = m(2c), 
zz m ±i — Z3 m ± (2c/3 + u), [so the puckering displacement from a layer is c/6 + u\ 
here m is any integer.] The total energy (per Al) is 

1 +1 

E= - [U(z m ) + V(z m+1 -z m )], (2) 

m— — 1 

where V(z' — z") is the Al interaction along the channel. [Note a small abuse in this 
approximation: since the channels are not straight, the actual Al-Al distance is not 
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just a function of z ©. Inserting z m into equation and Taylor expanding in u, we 
obtain 



Here E = U a + V(2c/3). Also, U\ = U'(2c/3), etc; inserting their values from 
equation 0, minimizing with respect to u, and using U\" <C Vi", we find the energy 
minimumn 



It can be seen this solution is favourable when a > 0; that is, the local mirror 
layer of the Al's in each channel lies in the TM-rich layer [as observed in real 
Al-TM structures.] Then, when two (equi)distant TM columns affect the channel, 
a — > a% + 02 in equation ® and the lowest E m - m is obtained when \<j\ + 02 =2. 
Thus the relaxed total energy of channel Al favours a 'ferromagnetic' alignment of 
orientations, as was claimed. 

Each channel breaks z-translation symmetry in that the mirror-layer Al sits at 
either z = or z = c, but not both. These channels interact, in a pattern 1 6 1 whereby 
four channels get occupied around a single TM column, as shown schematically in 
figure |3J;; viewed in projection, the mirror-layer Al atoms involved in the channels 
form a sort of bent cross. 

(Note added. The puckering instability is presumably related to some very low 
frequency phonon modes 1 15 1; the channels are also the conduit for the strong Al 
diffusion I16I . - ) 

5. Conclusion 

To summarize, the methods developed by Mihalkovic et al in a study of 
the 'basic' Ni-rich decagonal AI70C010N120 quasicrystal has successfully led to a 
structural understanding (at T = 0) of 'basic' Co-rich decagonal Al7oCo2oNiio, 
which is described by a completely different tiling, despite a very similar nearest- 
neighbor order. In this case, the key motif is a 12. 6A diameter decagon having a 
pentagonal core. 

We applied our approach to a unit cell the same size as W(Al-Co-Ni), a 
crystalline approximant of 'basic Co'. Many details of our prediction (most 
importantly the 13A Decagons, and the locations of their centres) are in agreement 
with the experimentally solved structure 1171 . However, W( Al-Co-Ni) has a few 
highly puckered sites which can not be captured in our approach which begins using 
a fixed-site list and relaxes afterwards. [Mihalkovic and Widom 1 10 1 have given a 
more thorough ab-initio investigation of W( Al-Co-Ni).] 

In general, the simulations with a fixed-site list - followed blindly - would lead 
to an imperfect description of Al atoms in ring 2.5 and ring 3 of the 13A Decagon. 
In simulations with relaxations (and in reality), many of these atoms displace from 
the layers, in a fashion which we have explained on the basis of the pair potentials. 




(3) 




(4) 
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